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We present deep Hubble Space Telescope ACS images of the post-starburst 
shell galaxy AM 0139-655. We find evidence for the presence of three distinct 
globular cluster subpopulations associated with this galaxy: a centrally concen- 
ts) ' trated young population (~ 0.4 Gyr), an intermediate age population (~ 1 Gyr) 
and an old, metal-poor population similar to that seen around normal galaxies. 
The g — I color distribution of the clusters is bimodal with peaks at 0.85 and 
1.35. The redder peak at g — I = 1.35 is consistent with the predicted color for 
an old metal-poor population. The clusters associated with the peak at g — I = 
0.85 are centrally concentrated and interpreted as a younger and more metal- 
rich population. We suggest that these clusters have an age of ~ 0.4 Gyr and 
solar metallicity based on a comparison with population synthesis models. The 
luminosity function of these "blue" clusters is well represented by a power law, 
<p(L)dL oc L~ 18 dL. Interestingly, the brightest shell associated with the galaxy 

■ harbors some of the youngest clusters observed. This seems to indicate that 
f--. ■ the same merger event was responsible for the formation of both the shells and 

the young clusters. The red part of the color distribution contains several very 
bright clusters, which are not expected for an old, metal-poor population. Fur- 

■ thermore, the luminosity function of the "red" GCs cannot be fit well by either 
^ ■ a single gaussian or a single power law. A composite (gaussian + power law) 

fit to the luminosity function of the red clusters yields both a low rms and very 
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plausible properties for an old population (with gaussian distribution) plus an 
intermediate-age population (with power-law distribution) of GCs. Hence, we 
suggest that the red clusters in AM 0139-655 consist of two distinct GC sub- 
populations, one being an old, metal-poor population as seen in normal galaxies 
and one having formed during a recent dissipative galaxy merger (likely the same 
event that formed the ~ 0.4 Gyr old clusters). 

Subject headings: galaxies: elliptical and lenticular — galaxies: clusters: general 
- galaxies: individual (AM 0139-655) — galaxies: interactions — galaxies: star 
clusters 



Introduction 



Globular clusters (GCs) are very useful probes of the dynamical and chemical assembly 
histories of galaxies. Many globular cluster systems in "normal" giant elliptical galaxies show 
a bimodal color distribution, providin g clear evidence for the presence of a "second event" 
during the formation of these systems (lAshman &: ZepJll998. and references therein! ) . Spec- 
troscopic studie s have revealed that both the blue and red GC populations are typ ically 



old (> 8 Gyr, Forbes et all l200ll : ICohen. Blakeslee. fc Cotel I2OO3I : IPuzia et al.ll2005h ) im- 



plying tha t the bimodality is ma inly due to a difference in their metallicity. The merger 
scenario of lAshman &: ZepJ (119921 ) predicted the presence of this bimodal distribution. How- 
ever, different theories were sug gested by variou s authors to explain th e nature of the sec- 
ond event. The merger model (jSchweizerl 1 19871 : lAshman fc Zepj 119921 ). predicts that the 
metal-rich cluster s are formed during major merger s of gas-rich galaxies. The "multi-phase 
collapse" model (IForbes. Brodie. fc Grillmairl 119971 ) implies that they form during a sec- 
ondary ^_cj3Uapj3e_pl^^ The accretion sce- 
nario (IForte. Martinez, fc Muzziol Il982l ; ICote. Marzke. &: Westl Il998l ) proposes that every 
galaxy is born with a GC system t hat has a median color ( metallicity) according to the 
color-magnitude relation of galaxies (le.g. Terlevich et al.lll999l ). The bimodality then arises 
due to the accretion of low-mass galaxies with metal-poor GCs by a larger galaxy with its 
metal-rich GCs. 



Intermediate - age GCs have been found in several galax ies (IGoudfrooij et al.ll200ll . 12007 



Puzia et al. 2002: Kundu et al. 2005: Chandar et al. 



20061 ). A few of these galaxies show 



features like shells or a disturbed morphology, indicating a recent merger e vent. Shells are 
sharp ly defined, arc-like features, mostly found around elliptical galaxies (IMalin fc Carter 
19831 ). though a few (Arp 215, Arp 227) have been found around disk galaxies ( Arplll9661 ). 
Two classes of models have been proposed to account for the shells. The internal models 
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( IFabian et al.lll980l ; IWilliams fe Christiansenlll985l ) suggest that shell stars form because of 
star formation within a giant shell of shocked interstellar gas. On the other hand, the interac- 
tion models suggest an external origin for the shell structure and interpret it as arising due to 
a merger. When an elliptical galaxy accretes a small companion galaxy, some vest iges of the 
event survive for a Gyr or longer after the merger in the form of ripples and shells (ISchweizer 
19801 ; lQuinnlll984j ). The morphology and alignment of the resulting shell system depends 
on the form of the potential of t he primary galaxy, the orbital parameters of the collisio n 



on tne iorm 01 tne potential 01 t ne primary galaxy, tne orpitai parameters 01 tne collisio n 
and the nature of the companion ( Hernquist fe Quinn 1988^ 19891 Dupraz fe Combesl 19861) . 
Typi cal dynamical ages of shells are 0.5 - 2 xlO 9 years ( Nulsen 19891 ; Hernquist fe Quinn 
19871). These s h ells, m ade up of sta rs, have been frequently observed around ellipticals. 
Malin fe Carter! Jl983h find 17% and ISchweizer fe Fordl Jl985l ) find 44% of field ellipticals 
to be surrounded by shells. The presence of shells thus seems to be a common outcome of 
interactions of galaxies. 

To test the merger scenario for GC formation, it is important to study galaxies that 
show obvious signs of a merger event in the form of shells, ripples, and some signature of 
having hosted a strong burst of star formation. One way of attempting this is to study 
E+A galaxies featuring shells. E+A galaxies were originally defined by the combination 
of an early-type galaxy morphology and the presence of a so-called "post-starburst" spec- 
trum containing strong Ba lmer absorption lines ( indicative of A-type stars) without strong 
[O n] or Ha emission lines ( Dressier fe Gunn 1983 1^1. The existence of strong Balmer absorp- 
tion lines indicates that these galaxies have experienced star formation around ~ 0.3 — 2 
Gyr ago while the absence of emission lines indicates that star formation must have ceased 
shortly thereafter. A large fraction (75%) of nearby E+A galaxies are located in the field 
(jZabludoff et al.lll996l ). Several of these show tidal features indicative of galaxy-galaxy in- 
teractions and mergers as the probable mechanism responsible for their formation. 



Long-slit spectroscopy of E+A galaxies ([Norton et al.ll200ll ) revealed that the young 
stellar populations are typically more central ly concentrated tha n the old population, but 
are not always confined to the inner 1 kpc. ISchiminovich et al.l (119941 ) discovered atomic 
hydrogen associated with the shells of NGC 5128. Subsequently, CO molecules were found 
in two shells surrounding this elliptical providing evidence that the molecular gas in merg- 
ers m ay not only fall into the n ucleus but can also be spread out to large galactocentric 



radii (ICharmandaris et al.ll2000i ). ICarter et al.l ( 



formation in 20 out of 100 shell galaxies in the 



988) found evidence for recent nuclear star 



Malin fe Carterl (119831 ) catalog, suggesting 



1 Subsequent studies revealed the presence of a disk component in many such galaxies, and many papers 



now use the term k+a, where the k stands for the spectral type of an old stellar population (e.g., 
1993h . 
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that the process of shell formation may involve some accretion or merger of gas-rich material 
which would account for the nuclear star formation. 

Numerical sim ulations of the sink i ng of a gas-rich satellite disk galaxy into a large 



elliptical galaxy by iKojima fc Noguchil (119971 ) show that the stellar particles from the disk 



stars of the satellite start making a clear shell structure after the first collision with the 
elliptical in case of a radial or a retrograde merger. This shell system lasts for ~ 1 Gyr or 
more. They also find an abrupt reduction in the star formation rate around the time of first 
passage of the satellite through the center of the elliptical due to the sudden decrease in gas 
density caused by the gravitational scattering by the elliptical potential. 

The study of GCs in such galaxies should yield important information regarding the 
nature of the E+A phenomenon. If numerous intermediate-age GCs with ages consistent 
with those of the post-starburst population in the nucleus were to be found, we could be fairly 
confident that the A stars were formed in a relatively vigorous and distinct star formation 
event, likely associated with a dissipative galaxy merger. Conversely, an absence of such 
GCs would indicate that no strong star formation was associated with the event that caused 
star formation to end ~0.5 to 2 Gyr ago. The latter would be consistent with a scenario 
where gas is removed from the galaxy, e.g. by ram pressure stripping. This paper focuses on 
a study of the globular cluster system of the early-type shell galaxy AM 0139-655. 



AM 0139-655 



AM 0139-655 (PGC 6240) is a giant e arly-type galaxy with bright shells that are irreg- 
ularly distributed around the main body (IMalin fc Carterlll983l ). The nuclear spectrum of 
AM 0139-655 shows evidence for a post-starburst (E+A) spectrum with very strong Balmer 
absorption lines (suggesti ng the presence of a substantial fraction of A-type stars) but no 
significant [On] emission (ICarter et al.lll988l ). AM 0139-655 has a large B.5 equivalent width 
(13.6 ± 1.0 A) and an active nucleus. Since AM 0139-655 shows both shell structure (in- 
dicative of a merger) and a post-starburst spectrum (signaling a burst of star formation 
and hence possible globular cluster formation in the recent past), it is a good candidate to 
explore whether any relations exist between its globular cluster system, its morphology and 
the merger event that gave rise to the shells. In this paper, we present a detailed photomet- 
ric study of the GCs of AM 0139-655. The general properties of AM 0139-655 are listed in 
Tabled 
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3. Observations and Data Reduction 

AM 0139-655 was observed with HST on July 16 and 17, 2004, as part of General 
Observer program 10227, using the Wide Field Channel (WFC) of ACS. It was observed 
for a total of 6 orbits in the continuous viewing zone (CVZ) of HST with total integration 
times of 21440 s in the F475W and 7300 s in the F814W filter. The observations consisted 
of several long exposures at different dither positions as well as a few short ones to prevent 
saturation of the galaxy center. The data were processed with the ACS on-the-fly pipeline 
which included dark and bias subtraction and flat fielding. Each flat-fielded image was 
carefully checked for satellite trails and saturated pixels in the central region of the galaxy. 
These pixels were flagged and masked out. The individual images in ea ch band were co- 
added using the PyRAF task MULTIDRIZZLE feoekemoer et alJl2002h . This procedure 



removed all the cosmic rays and corrected for geometric distortion. The combination of the 
short and long exposure images enabled us to get good resulting multidrizzled images with 
unsaturated centers. The final combined F475W image is shown in Fig. [TJ 



4. Photometry 



The search for sources was performed on an image made from the summed F475W+F814W 
image in the following manner. The F475W and the F814W images from multidrizzle were 
added together. To improve the detection and the accuracy of the photometry, especially for 
the faint sources embedded within the bright galaxy background, the underlying galaxy light 
was first modeled. Extensive experimentation showed that the application of a circular me- 
dian filter gave the best resu lts. Ob ject de tection was done by using the DAOFIND routine 
in the DAOPHOT package (jStetsonl 119871 ) on an image obtained by dividing the combined 
image by the square root of the median filtered image (thus providing uniform shot noise 
characteristics over the whole frame). The advantage in using the g+I image rather than the 
individual images for source detection is the greater depth reached by the combined image. 
Also, the photometric zero-point of g+I is significantly less dependent on the color of the 
source than the individual g or / images, which helps during the completeness correction 
(see Goudfrooij et al. 2007 for a detailed discussion). The detection threshold was set at 
4cr abo ye the background. Typical half-light radii of globular clusters are in the ran ge 1- 
20 pc fcundu k Whitmorej l200ll : Ivan den Bergh k Mackevl 12004! ; Ijordan et al-lfeood ). At 
the distance of AM 0139-655, the scale is 25 pc per ACS/WFC pixel. Thus, we expect the 
globular clusters to be nearly unresolved point sources on the image. 

Smooth, median-filtered model images were created in each of the two bands from the 
corresponding multidrizzle combined images in a manner similar to the one described above. 
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The model was subtracted from the corresponding original images to get a residual image 
in which the galaxy light has been removed to a great extent and the compact objects are 
clearly visible. To confirm that the subtraction procedure did not alter the error statistics of 
the photometry, we compared the magnitudes obtained from the subtracted and the original 
images for a few bright sources in the frame. The differences between the two methods 
were well below the photometric errors. Aperture photometry was carried out on these 
residual images through an aperture of 3 pixel radius (0'/15). In order to discard obvious 
spurious detections (e.g. edges of the image), objects with magnitude errors greater than 
0.4 mag were rejected at this point. Aperture corrections from 3 pixel (0'.'15) to 10 pixel 
(C/5) radius were determined using a few bright poin t sourc es in each band. The corrections 
from 0".5 to infinity were taken from ISirianni et al.l (120051 ). The total aperture corrections 
were 0.470±0.009 and 0.517±0.010 in F475W and F814W respectively. The transformation 
from the instrumental magnitudes (STMAG) to SDSS g- and Cousins /-band was done using 
the SYNPHOT package in STSDA S . Synt hetic spectra of stellar types ranging from BOV to 
K3III from the iBruzual fc Charlotl (119931 ) library were convolved with the SDSS g and the 
Cousins / filter to obtain magnitudes relative to Vega, and with the F475W and the F814W 
filters to yield magnitudes in the STMAG system. Fitting polynomials to get the relation 
between the two systems resulted in the following transformations: 



g = FA75W + (0.4530 ± 0.0006) 

+ (0.039 ± 0.001)(F475iy - F8UW) 

+ (0.0010 ± 0.0003) (FA7hW - F8UW) 2 

-(0.0040 ± 0.0003) (FA75W - F8UW) 3 
I = F8UW - (1.258 ± 0.002) 

+ (0.015 ± 0.002)(F475W - F8UW) 

-(0.002 ± 0.003)(F475iy - F8UW) 2 

-(0.001 ± 0.001) (F475iy-F814W0 3 (1) 



4.1. Globular Cluster Candidate Selection 

The photometry lists with the magnitudes in g and / obtained in the preceding section 
also included foreground stars and extended objects, in addition to cluster candidates. To 
remove these contaminating sources we applied two selection criteria based on the FWHM 
and the compactness of each source. We obtained the direct FWHM of each of the sources 
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in the photometry list in each band using the task IMEXAM within IRAF@ in batch mode. 
Since the FWHM of bright stars in our frame is ~2.0 pixels and we expect all the cluster 
candidates to be nearly unresolved point sources, we discarded all sources that did not satisfy 
the criterion 1.5 < FWHM < 3.0 pixels. We also added a compactness criterion wherein only 
sources that had the difference in the magnitudes between aperture radii of 2 pixels and 5 
pixels lie between 0.2 mag and 0.6 mag were retained. This ensured (to a great extent) that 
the contamination from background galaxies was minimized. To study the radial extent and 
distribution of the GC candidates so obtained, they were binned into annuli 400 pixels wide, 
starting from the center of the galaxy out to the corners of the ACS frame. An effective 
area was calculated for each annulus namely the total area of the pixels that lie within the 
annulus which had a non-zero value of intensity. This ensured that the masked-out regions 
like the chip gap were not considered while computing the area since these areas were also 
excluded during cluster finding. We plot the logarithm of the number of sources per kpc 2 as 
a function of the projected distance from the center of the galaxy and find that it reaches an 
asymptotic value at a projected radius of about 40 kpc (see Fig. [2]). Given the large distance 
of the galaxy and hence, the small angular size of AM 0139-655, it was possible to use the 
outer parts of the image to estimate the contamination due to compact background sources. 



From the lBurstein & Heilesl (119 82) model, we find the Milky- Way foreground extinction 
Ay = 0.00. The resulting g versus g—I color-magnitude diagram for the candidate globular 
clusters is shown in Fig. [3l The rectangle delineates the sources that meet our criteria 
for globular cluster candidates. Using the GALEV simple stellar population (SSP) models 



(lAnders &: Fritze-v. Alvenslebenll2003l ). we determine the color range of < g — / < 2.0 to 
include all clusters over the range of ages from 10 7 to 10 10 years and metallicity z values 
from 0.0004 to 0.05 (where z = 0.02 represents solar metallicity). All sources with colors 
redder than 2.0 are expected to be late-type foreground stars or red background galaxies and 
are rejected in the subsequent analysis. The final list, thus obtained, contains a total of 217 
globular cluster candidates. Typical photometric errors in g were 0.01, 0.03, 0.15 and >0.2 
mag for magnitude intervals 22-24 mag, 24-26 mag, 26-28.75 mag and fainter than 28.75 
mag respectively. The brightest fifty GC candidates are listed in Table [2j 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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4.2. Completeness Corrections 

Completeness tests were performed on the combined g+I image, i.e., the one used for 
detection of cluster candidates. Artificial objects were added in batches of 100 for five 
background levels and several magnitude intervals. The background levels were computed 
from the g+I image in the following manner. The counts in the outermost parts of the image, 
free from any galaxy light contamination were taken as Background level #1. The counts 
in the innermost region of the galaxy where clusters are detected were taken as Background 
level #5. This difference between the faintest and the brightest background levels was then 
divided into five equal intervals in log space to get the five mean background levels. The 
radial intensity profile of the artificial input globular clusters was determined by fitting PSFs 
to the cluster candidates in the original image. The image was divided by the square root of 
a median-filtered version of itself and source detection was performed on the resultant image, 
with the same procedure as used for the initial source finding. Further, the ranges of the 
FWHM, compactness and maximum permissible error in the derived magnitudes were also 
selected in exactly the same manner as that for the actual source photometry. The resultant 
completeness curves are shown in Fig. HI 



5. g—I Color Bimodality 

To analyze the color distribution of globular clusters in AM 0139-655, we start by 
computing the non-parametric Epanechniko y-kernel probabil ity density function Pobj of all 



objects inside the analysis radius of 40 kpc (jSilvermanl Il986[ ) . Given the photometric com- 
pleteness of our data "P bj is the closest approximation to the real color distribution of globular 
clusters. The width of the adaptive kernel linearly depends on the statistical color uncer- 
tainty of every individual object, so that the initial noise statistics of the color measurements 
is translated into the final probability density estimate. To account for uncertainties due to 
the limited sample size, we bootstrapped this process 100 times and derive 90% confidence 
limits of the resulting probability density function. 

With the same technique we also computed the probability density function of the 
background field Vbkg derived from the detections. We then scaled the background density 
estimate to the object function V h] by area. We note that the error associated with the 
scaling of the function is included in our bootstrap routine and propagated to the uncertainty 
of the final density estimate. We finally derive the color distribution of globular cluster 
candidates by statistical subtraction: 



PGCC = ^obj - ^bkg 



(2) 
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where A is the ratio of areas encompassed by the object versus background regions. Fig- 
ure [5] shows the background corrected function. The final color distribution appears clearly 
extended and, taking the 90% error margins into account, bimodal with peaks at g — I = 
0.85 and g — I — 1.35. A rather ext ensive red wing to the red peak is present, whic h will 
be discussed in Sect. 7. KMM tests ( jMcLachlan fc Basfordlll988l ; lAshman et al.lll994l ) were 
also run on the background subtracted g—I color distribution. The KMM test for two com- 
ponents resulted in peak values consistent with those found using the Epanechnikov kernel. 
The test for three components resulted in a third peak with g — I = 1.45 but without any 
clusters assigned to it. 

Although we do not know the effects of possible internal extinction on individual GCs, 
the fact that the blue peak of the GC color distribution found in AM 0139-655 at g—I = 1.35 
is consistent with V — I = 0.95 (for 14 Gyr old, metal-poor clusters from SSP models) which 
agrees with the blue peak of GCs in the Milky- Way indicates that on average the internal 
visual extinction is very small. At the same time, this also support s our assumption that the 
Milky- Way foreground extinction is near zero, as predicted by the iBurstein &: Heilesl (Il982l ) 
model. 



5.1. Spatial Distribution of the Blue and Red Subpopulations 

As shown in the preceding section, the color distribution of the clusters has a bimodal 
form with peaks near g — I = 0.85 and 1.35. The spatial distributions of the two subpopu- 
lations were derived by dividing the clusters into two subpopulations based on their colors. 
Clusters which lie in the range < g — I< 0.95 were designated as the blue population and 
those within 1.15 < g — I< 2.0 were designated as the red population. The region 0.95 < 
g—I < 1.15 was avoided to minimize mixing of the two populations. This yielded 84 blue 
clusters, 114 red clusters and 19 clusters that lie in the region of mixing. 

Figure [2] shows the radial distribution of the blue and red samples. The bins were chosen 
in such a way that there were at least five clusters in each bin in the region within 40 kpc. 
It is clear from the figure that the blue clusters are more concentrated towards the center 
and their surface density falls as we go away from the center. On the other hand, the surface 
density of the red clusters is somewhat low in the center and may peak between 5-10 kpc 
from the center of the galaxy. The surface density distribution of the 19 GCs that lie in the 
intermediate color range is also shown in Fig. [2j This distribution appears to peak at the 
center. Adding all these clusters to the blue clusters would not change the trend for the blue 
clusters while adding them to the red clusters would tend to somewhat dilute the effect in 
the central regions. 
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We used the SSP models described in lAnders fc Fritze-v. Alvenslebenl ( 120031 ) to esti- 
mate the ages of the two populations. These models include continuum and line emission 
and take into account the important contribution of gaseous emission to broad-band fluxes 
and their strong metallicity dependence during the very early evolutionary stages of star 
clusters. The top panel of Fig. [6] shows the dependence of g — I color on the age of the 
cluster. From the figure, we see that the red peak in the color distribution at g — I = 1.35 
can be explained by clusters of age 14 Gyr and metallicity Z= 0.0004 which agrees very well 
with the colors expected for the old, metal-poor globular cluster population that forms the 
blue peak in normal ellipticals or spirals. The blue peak of the color distribution at g — I = 
0.85 can be explained by a younger population of clusters of solar metallicity with an age ~ 
0.4 Gyr. Note that this age estimate is fairly independent of the metallicity: See Fig. [6j This 
is consistent with the E+A nature of AM 0139-655. We assume solar metallicity for these 
clusters since several spectroscopic studies of luminous intermediate-age globular clusters in 
a number of merger remna nts have sho wn that they have near -solar metalli cities (NGC 7252 



(ISchweizer & Seitzeilll998h: NGC 3921 (Schweizer et alJl2004h : NGC 1316 (IGoudfrooii et al. 



20011 ) and NGC 5128 (IPeng et al.l 120041 )) . HST imaging studies of the same merger rem 
nants show that these youn g clusters are more concentrated towards the galaxy center than 
the old metal-po or clusters ( Miller et al. 19971; Schweizer et al. 1996 : Goudfrooij et al. 2001 



Peng et al.l 120041 ) , as predicted by lAshman &: Zepj (jl992h who modeled the formation of 
clusters in a dissipative galaxy merger. 

Figure [7J shows the projected spatial distribution of the blue and the red clusters while 
Fig. [S] shows the g—I color map obtained by dividing the g image by the / image. In Fig. [HI 
note the presence of a blue ring-like structure surrounding the center and the blue clusters 
present along the ring. Another interesting fact is that the brightest and the most distinct 
shell (lying between 14.8-16.7 kpc to the East of the galaxy center) harbors seven blue 
clusters (Fig. [2]), all of which fall within the color range 0.70 < g—I < 0.91, consistent with 
SSP model ages between 0.25-0.4 Gyr. Thus, it appears that all the clusters associated with 
this shell were formed during the latest burst of star formation. We suggest that the same 
event that formed the shell structure also gave rise to the young globular cluster population. 



6. Luminosity Function of AM 0139-655 

To derive the luminosity functions (LFs) of the clusters, we proceeded in the following 
manner. A completeness correction was applied to each cluster depending on its magnitude 
and background level. This was determined from the curves shown in Fig. 0] using bilinear 
interpolation. The clusters were then divided into two groups depending on their distance 
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from the galaxy center. The clusters within a projected distance of 40 kpc from the center 
of the galaxy were considered as actual GC candidates and the ones farther were considered 
as contaminants (foreground stars and compact background galaxies). The LFs of the GC 
candidates were then corrected for background contamination using the smoothed and scaled 
LFs of these objects. The LFs obtained after completeness and background contamination 
correction and the uncorrected LFs in the /-band are shown in Fig. [TUl 

Given the distance modulus of 35.3 for AM 0139-655, the globular clu ster turn-ove r 



magnitude for old, metal-poor GCs is expected at / = 27.2 (for Mj = -8.1) (IHarria Il996l ) . 
As seen from the figure, the data does not quite reach the peak of the turnover. We analyze 
the luminosity functions of the blue and red clusters separately. 



6.1. Luminosity Function of the Blue Clusters 



As mentioned above in Section 5.1, all clusters with colors in the range 0.0 < g— / < 0.95 
are considered as blue clusters. Figure [10] (middle panel) shows the luminosity function of 
these blue clusters. A power-law fit to the LF of the blue clusters yields </>(L)dL oc L~ a dL 
with a = —1.8 ± 0.1, consistent with the results of previous LF studies of young clusters in 
mergers and merger remnants which find power- law indices in the range —2.1 < a < — 1 . 7 
dwhitmore & Schweizerlll995l : iMeurer et aDll995l : Iwhitmore et alill999l : Iwhitmorel 120031 ) . 



6.2. Luminosity Function of the Red Clusters 

Clusters within the color range 1.15 < g — I < 2.0 were classified as red clusters and 
their luminosity functions derived (Fig. [TO], bottom panel). Attempts to fit their luminosity 
function by either a single gaussian or a single power law yielded large residuals. Upon close 
examination of the color-magnitude diagram, we found a significant number of very luminous 
red clusters. If the red clusters were attributed solely to an old, metal-poor population, with 
a mass function similar to that in our Galaxy, the brightest clusters would not be expected 
to be brighter than M g = —10.0 (g = 25.0), yet we find an excess of clusters as bright as g 
= 23.5. In addition, the color distribution of the red clusters shows a broad tail towards red 
colors (Fig. [5]) which suggests the presence of more than one sub-population. 
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6.2.1. Modeling the Composite Luminosity Function 



Spectrosco pic studies of inter mediate- age clusters in NGC 1316 (IGoudfrooij et al.ll200ll ) 
and NGC 5128 (IPeng et al.ll2004l ) have shown that they have near-solar metalli c ities. Simu- 
lations of the color- magnitude diagrams of evolving clusters by I Whitmore et al.l (119971 ) show 
that solar metallicity clusters in the age range of 1-2 Gyr are the most difficult to identify 
since they have nearly the same colors as the old, metal-poor GC subpopu lation. However, 



they a re significantly (2-2.5 mag) brighter than the old ones for a given mass. IWhitmore et al. 
(119971 ) also find that solar metallicity clusters older than ~ 2 Gyr will be redder than the 



old, metal-poor cluster population. Using the SSP models shown in Fig. El we see that the 
range of colors corresponding to the red cluster distribution in AM 0139-655 can indeed be 
explained by a metal-poor population about 14 Gyr old, but also by a metal-rich younger 
population about 1-1.5 Gyr old. To test whether our observed distribution may indeed be 
the result of two distinct populations of the same color but different ages and metallicities, 
we model the observed luminosity function as a combination of two populations; one old, 
metal-poor population which is consistent with a gaussian magnitude distribution and the 
othe r an intermediate-age popul ation of solar metallicity, which can be described by a power 



law (je.g. Goudfrooij et al. lion ). 



To constrain the possible relative contributions of the two populations, we consider the 
specific frequency of GCs, 

S N = N GC W 0A{Mv+15) (3) 

i.e., the number of star clusters per gala xy luminosity normalized to an absolute F magnitude 
of —15 ([Harris fc van den Berghl 119811 ). To estimate the contribution of the old, metal-poor 
cluster population to the mixed luminosity function, we need to first consider the possible 
Sn values of the progenitors and their V magnitude. Assuming AM 0139-655 was formed by 
a merger of two galaxies, the progenitor galaxies could have been two spirals or one elliptical 
and one spiral. Since there is no a priori obvious way of determining the nature of the 
progenitors, we consider various scenarios. 

The specific frequency is known to increase systematically along the Hu bble sequence , 
from 0.5 ± 0.2 for Sc spirals to 2.6 ± 0.5 for ellipticals outside of clusters (jHarrisI Il99ll ). 
Hence, we consider four values of Sn (0.25, 0.5, 1.0 and 3.0) covering the range from spirals 
to ellipticals to describe the possible specific frequencies of the progenitors. We use the V 
magnitude of AM 0139-655 to estimate the combined V magnitude of the progenitors. This 
will be a rough estimate and will tend to overestimate the luminosity since star formation 
during the merger would brighten the merger remnant. However, the V magnitude of AM 
0139-655 is not listed in the literature. Hence, we first calculate the total g and / magnitudes 
of AM 0139-655 from our images using ellipse fits and calculate g—I from these magnitudes. 



13 



We get g = 13.64; 1= 12.34 and hence, g-I = 1.30. 

Using this, we find the luminosity-weighted age of 1 Gyr for a solar metallicity population 
at this g—I from Fig. [6J A 1 Gyr solar metallicity population has V— I = 0.89 and hence V 
= 13.23 according to the GALEV models. We then compute the fading in the Uband from 
this age to an age of 14 Gyr from the lower panel of Fig. [6j We use luminosity-weighted age 
at solar metallicity as we are attempting to separate the intermediate-age clusters (which 
should have about solar metallicity) from the old, metal-poor GC population. Though a 
SSP is not a very good approximation to a galaxy, it is sufficient for our purpose since 
we are interested in getting a rough estimate of the fading of the galaxy. Using the faded 
value of V — 15.84 (My = -19.45), in the equation for Sn, we calculate the number of 
old globular clusters expected for the four values given above. The estimated ga ussians 



for the old, metal-poor population we re calculated using Mj = —8.1 and a — 1.2 ([Harris 



19961 ; iBarmby. Huchra. fc Brodidl200ll ) for each Sn- Fig. [TT] shows the observed luminosity 
function of the red clusters and the estimated gaussians for the four values of Sn- Subtracting 
the contributions to the gaussians in each case and examining the residuals, we find that the 
residuals in each case can be fit by power laws with power-law indices a between -1.6 and 
-1.85 consistent with values found for intermediate-age cluster populations. The rms values 
of the difference between the power law fits and the histograms in the right panel of Fig. [11] 
are 4.08, 2.88, 3.28 and 3.02 for Sn = 3.0, 1.0, 0.5 and 0.25 respectively. The above analysis 
suggests that the red cluster population in AM 0139-655 is in fact a combination of two 
populations of different ages. To conclude, we have identified three populations of GCs in 
AM 0139-655: one blue, metal-rich population of young clusters and a red population that 
can be split into two populations, one old (14 Gyr) metal-poor and one of intermediate age 
(1-1.5 Gyr) and solar metallicity. Since the two populations contributing to the red peak had 
similar colors, the KMM test could not identify them as separate populations. The estimates 
of ages, metallicities and the number of clusters (corrected for background contamination 
and completeness) belonging to the three sub-populations are given in Table [3] for the various 
possible scenarios. However, spectroscopic observations of the brightest candidates from the 
three populations will be necessary to derive exact metallicities and verify our tentative 
interpretation. 



6.3. Masses of the Clusters 

It is illustrative to estimate the masses of the bright clusters formed during the merger 
and compare them to the masses of the most massive clusters in our Galaxy (u Cen with a 



mass 5 x 1O 6 M ) and M31 (Mayall II, Gl), with a mass (1.4 - 1.7) x 1O 7 M flMevlan et al. 
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200l|). 



Fig. [TJ] shows curves of constant mass superimposed on the color-magnitude diagram. 
To derive these masses, a grid was constructed for the observed magnitude and color range. 
These magnitudes were then faded to an age of 14 Gyr. Masses were a ssigned to each cluste r 
by assuming the value of M/Ly to be 4.1 (at 14 Gyr), as found by iMeylan et al.l (119951 ). 
We find that several clusters in AM 0139-655 will have masses greater than that of the most 
massive cluster in our galaxy and comparable to the most massive cluster Gl in M31 when 
they age to 14 Gyr without subsequent mass loss due to stellar e volution since this mass 
loss is small beyond 1 Gyr (lAnders &: Fritze-v. Alvenslebenl 120031 ). The clusters will lose 
stars during their dynamical evolution. Their stellar-dynamical mass loss, however, cannot 
be evaluated with any confidence without orbital information. 

The masses derived in this way have two cave ats. First, the mas s-to- light ratio is 
assumed to be equal to that of u Cen. Various studies (IMeylan et al.lll991l) find a wi de range 
(1.2 to 5.2) of mass- to-light ratios fo r GCs. According to iBruzual fc Charlotl (120031 ) models 
computed using the IChabrierl (120031 ) IMF and lower and upper mass cut-offs mi = 0.1M Q 
and m;/ — 100M Q , M/Ly of a 14 Gyr solar metallicity population is 4.3 while that of a 
1/50 Z Q population of the same age is 2.0. Hence, a mass-to-light ratio of 4.1 is a good 
approximation for solar metallicity clusters. In the absence of any spectroscopic data on 
these clusters, these values can be taken as rough estimates. Secondly, the effects of internal 
reddening (though expected to be small from Sect. 5), have not been taken into account. 
This is not expected to change the mass ranges much since dereddening would tend to shift 
the points to brighter magnitudes and bluer colors (upward and to the left in Fig. (T2J), which 
is nearly parallel to the curves of equal mass. 



7. Nature of the Progenitors 



In a recent study of GC systems around six shell ellipticals, ISikkema et al.l (120061 ) typi- 
cally find bimodal distributions similar to those of normal ellipticals. They find some indica- 
tions of the presence of a younger intermediate-age globular cluster population in two of their 
galaxies, but the number of these intermediate-age clusters is small as compared to the old 
clusters. However, in the case of AM 0139-655, we find that the cluster system is dominated 
by the young and the intermediate-age clusters with a relatively small contribution from the 
old, metal-poor clusters. In order to constrain the number of old, metal-rich clusters that 
might have been inherited from the progenitors, we note that there is some evidence for a 
small peak which is consistent with a redder, metal-rich sub-population (at g—I = 1.8) as 
expected for a sub-solar metallicity from the SSP models. The small number of GCs in this 
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minor peak suggests the absence of a significant bulge in either precursor. Additionally, for 
the four values of Sjv considered for modeling the LF of the precursors, the rms values point 
out that at least one of the progenitors was not an e lliptical, but had small 5W values < 1, 



which is consistent with Hubble types of Sb or later (IChandar et al.ll2004l ; iGoudfrooij et al. 



20031 ) which have small bulges. Both the above findings seem to suggest that the galaxies 



involved in the merger that formed AM 0139-655 did not have a large bulge generally seen 
in early-type galaxies but were most likely late-type spirals with small bulges. A merger 
involving late-type gas-rich spirals could have formed one population of new GCs about 1 
Gyr ago in a starburst a t first pericenter and another one 0.4 Gyr ago at final merging 
flMihos fc Hernquistl[l99~6] ). 

From the fact that we see evidence for a young and intermediate-age component of 
clusters in AM 0139-655, we can safely assume that at least one of the galaxies involved 
in the recent interaction must have been rich in gas since there were two episodes of star 
formation within the time periods for shell survival. This is also supported by the fact that 
AM 0139-655 itself has strong Balmer absorption lines, which suggests that at least one of 
the progenitors must have been gas-rich and hence could not have been an early-type with 
a large bulge. 



8. Summary and Conclusions 

Deep ACS images of the post-starburst shell galaxy AM 0139-655 have been analyzed. 
We find a total of 217 GC candidates associated with this galaxy. The g— I color distribution 
of the clusters shows two distinct peaks, one at g — I = 0.85 and the other at g — I = 1.35. 
We compare the colors with SSP models and find that the red peak coincides with the 
color expected of the old, metal-poor population. The blue peak corresponds to a young 
population (0.4 Gyr) assuming solar metallicity. The luminosity function of the blue clusters 
can be fit by a power law of index a ~ -1.8±0.1, as is expected for young clusters. The red 
part of the distribution shows a long tail and the presence of several luminous clusters. We 
show that the luminosity function of the red clusters can be modeled as a combination of 
two subpopulations: An old population with a gaussian distribution and an intermediate-age 
population with a power law distribution. Thus, three subpopulations of globular clusters 
are clearly distinguishable in AM 0139-655. The number of clusters in the old, metal-rich 
subpopulation that is commonly associated with "normal" early-type galaxies is very small 
relative to the three aforementioned subpopulations. The combination of this dearth of old, 
metal-rich clusters and the fact that the LF of the old population is best fit by a small value 
of Sn supports the notion that any bulge in the progenitor galaxies must have been small. 
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The ages of the two younger GC subpopulations are comparable to typical lifetimes of 
shell structures formed during simulations of galaxy mergers, and the brightest shell actually 
seems to host several young GCs. This strongly suggests that both the shells and the young 
and intermediate-age GCs were formed during a dissipative merger. It thus seems fairly 
likely that the E+A signatures of this galaxy are due to this recent merger. These findings 
lend credence to the "merger scenario" for the formation of metal-rich GCs and their galaxy 
hosts. 
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Table 1. General properties of AM 0139-655. 



Parameter 


Value 


RA*(J2000) 


lh 41m 30.98s 


Dec*(J2000) 


-65° 36' 55.44" 


Hubble Type 


SO 


Vhel 


8216 ± 15 km s" 1 


v hel (corr) 


7936 km s" 1 


velocity dispersion 


249 ± 46 km s" 1 


Distance** 


105.8 Mpc 


m-M 


35.3 


scale 


1" = 512 pc 




1 ACS pixel = 25pc 


M B 


-20.57 


galactic latitude (b) 


-50.7° 




0.00 



Trom ACS images, this work. 
'*Using H = 75 km s _1 Mpc _1 



Burstein fc Heiled (119821 ) 



Note. - All other parameters are taken 
from LEDA Qhttp://leda.univ-lyonl.fr/D 
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Table 2. Positions and photometry of the 50 brightest globular cluster candidates 



RA ( 


;2ooo) 


Dec (2000) 


/ (mag) 


9-1 


(mag) 


r*(") 


r** (kpc) 


1 


41 


iq kk 

0U.00 


PK IP A K QQ 

-oo oo 4o. yy 


20.87 


± 


0.01 


117 
1.1 ( 


± 


0.02 


q ao 
y.oz 


K Q 
O.U 


1 


41 


iq qi 
OU.UO 


PK IK KA A1 
-00 00 04.40 


20.94 


± 


0.01 


1 1 1 
1.11 


± 


0.02 


fi1 OQ 

oi. zy 


1A A 
01.4 


1 


41 


ol.lo 


PK IP KA SO 
-00 00 01. oU 


21.12 


± 


0.01 


n Q7 
u.y < 


± 


0.02 


q 77 
O. 1 I 


1 Q 

i.y 


1 


41 


07 10 
Z 1 .OZ 


PK 17 A 1 07 
-00 ( 41. Z 1 


21.18 


± 


0.01 


1 1A 
1.04 


± 


0.02 


ki 1 n 

01. 1U 


OP 1 
ZO. 1 


1 


41 


1A 7Q 

or. 1 y 


PK IP KA QI 
-00 00 01. UO 


21.61 


± 


0.01 


n fi.fi 

U.oO 


± 


0.02 


p fia 

O.Oo 


1 A 
0.4 


1 


41 


11 KQ 

oo. oy 


pk ip on nn 
-00 00 zu.uu 


21.68 


± 


0.01 


1 on 
1 .zu 


± 


0.02 


1% QP 

oo.yo 


1 Q Q 

iy.y 


1 


41 


°.1 AO 

01 .4Z 


PK IP KK KQ 

-oo oo oo. oy 


21.75 


± 


0.03 


1 PO 
1 .OZ 


± 


0.07 


O 70 
Z. 1 z 


1 A 
1.4 


1 


41 


AT PQ 

41 .oy 


PK 17 1 K A K 
-00 ( 10.40 
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± 


0.01 


n pq 
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± 


0.02 


pq ir\ 
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IK K 
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1 


41 
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PK IP AP K1 
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± 
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± 
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lo. 14 


Q 1 
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1 


41 


iq io 
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PK IP K1 1 K 
-00 00 00.10 
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± 
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n so 
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± 
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A PQ 
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O A 
Z.4 


1 


41 


OA QP 

z4.yo 


pk ip on IK 
-00 00 ZU.00 


22.46 


± 
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n an 

U.oU 
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OP 
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Z4. iy 


PK 17 OQ 70 

-oo o ( zy . i u 
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0.04 


0.31 


± 


0.05 


3.45 


1.7 


1 


41 


30.16 


-65 36 54.38 


23.72 


± 


0.03 


0.96 


± 


0.04 


5.20 


2.6 


1 


41 


29.12 


-65 37 30.60 


23.81 


± 


0.02 


1.30 


± 


0.03 


36.99 


18.9 


1 


41 


30.15 


-65 36 59.70 


23.85 


± 


0.02 


1.40 


± 


0.04 


6.68 


3.4 


1 


41 


36.63 


-65 36 58.83 


23.95 


± 


0.02 


1.30 


± 


0.03 


35.15 


18.0 


1 


41 


36.02 


-65 36 46.84 


23.97 


± 


0.02 


0.83 


± 


0.03 


32.40 


16.6 
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Table 2 — Continued 



RA (2000) Dec (2000) /(mag) g - 1 (mag) r*(") r** (kpc) 



1 


41 


35.53 


-65 36 49.60 


24.02 


± 


0.02 


0.82 


± 


0.03 


28.78 


14.7 


1 


41 


26.32 


-65 36 59.17 


24.06 


± 


0.02 


1.37 


± 


0.03 


29.06 


14.9 


1 


41 


35.04 


-65 36 47.56 


24.10 


± 


0.02 


1.35 


± 


0.03 


26.35 


13.5 


1 


41 


21.85 


-65 36 2.19 


24.10 


± 


0.02 


1.23 


± 


0.03 


77.69 


39.8 


1 


41 


29.50 


-65 36 52.78 


C\ A 1 C\ 

24.12 


i 

± 


0.03 


0.85 


l 

± 


0.03 


9.51 


A O 

4.8 


1 


41 


22.49 


-65 37 31.88 


24.15 


± 


0.02 


1.97 


± 


0.04 


63.96 


32.7 


1 


41 


32.33 


-65 36 49.64 


24.16 


± 


0.02 


0.94 


± 


0.03 


10.18 


5.2 


1 


41 


29.73 


-65 36 48.36 


24.17 


± 


0.03 


1.39 


± 


0.04 


10.49 


5.4 


1 


41 


42.40 


-65 37 29.54 


24.26 


± 


0.02 


1.74 


± 


0.04 


78.53 


40.2 


1 


41 


31.52 


-65 37 0.55 


24.26 


± 


0.03 


0.93 


± 


0.04 


6.11 


3.1 


1 


41 


33.11 


-65 37 37.96 


24.35 


± 


0.03 


1.36 


± 


0.04 


44.52 


22.8 


1 


41 


24.02 


-65 36 42.34 


24.40 


± 


0.03 


1.35 


± 


0.04 


45.02 


23.0 


1 


41 


29.77 


-65 37 9.49 


24.42 


± 


0.03 


1.38 


± 


0.04 


15.90 


8.1 


1 


41 


29.56 


-65 37 1.17 


24.47 


± 


0.03 


0.85 


± 


0.04 


10.45 


5.3 


1 


41 


20.68 


-65 36 45.31 


24.50 


± 


0.03 


1.64 


± 


0.04 


64.55 


33.0 


1 


41 


28.81 


-65 37 3.23 


24.51 


± 


0.03 


0.73 


± 


0.04 


15.49 


7.9 


1 


41 


38.81 


-65 37 12.53 


24.56 


± 


0.03 


1.54 


± 


0.04 


51.41 


26.3 


1 


41 


34.14 


-65 35 44.84 


24.57 


± 


0.03 


1.20 


± 


0.04 


73.27 


37.5 


1 


41 


24.77 


-65 36 31.55 


24.57 


± 


0.03 


1.72 


± 


0.05 


45.29 


23.2 



*projected radius derived using the position of the center at RA = lh 41m 
30.98s and Dec = -65° 36' 55.44" 

"Using Hq = 75 km s^Mpc" 1 
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Table 3: Properties of the three sub-populations of GCs 





Age (Gyr) 


Z 


N 




Old 


14 


0.0004 


25 (for S N 
12 (for S N 
6 (for S N ; 


= 1.00) 
= 0.50) 
= 0.25) 


Intermediate 


1 1.5 


0.02 


81 (for S N 
93 (for S N 
99 (for S N 


= 1.00) 
= 0.50) 
= 0.25) 


Young 


0.4 


0.02 


50 






Fig. I.— ACS image of AM 0139-655 in F475W 
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Fig. 2. — Surface density of the globular cluster candidates as a function of the distance 
from the galaxy center. The black solid line shows all the clusters, the blue line shows the 
distribution of the blue clusters and the red line shows the distribution of the red clusters. 
The surface density of clusters lying in the color range 0.95 < g—I < 1.15 (between the blue 
and the red clusters) is shown by the green line. The bin width is 10 kpc except for the 
inner 20 kpc for the blue and the red clusters where it is 5 kpc. The blue clusters seem to 
be more centrally concentrated as compared to the red clusters. 
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Fig. 3. — g versus g—I color-magnitude diagram for AM 0139-655 globular cluster candidates. 
The open squares denote GC candidates present within the inner 40 kpc and the crosses 
denote those present beyond 40 kpc. The dotted lines delineates the region of GC selection. 
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(g + D 

Fig. 4. — Completeness correction curves for the combined g and I band for five background 
levels. From left to right: 2000, 980, 480, 200, 120 e~ per pixel for an effective (/-/-/exposure 
time of 1822 s. 
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Fig. 5. — GC color distribution for AM 0139-655. The open histogram in the top panel repre- 
sents the observed color distribution and the hatched histogram represents the distribution 
after subtraction of the background. The background distribution is shown in the lower 
panel. The solid red line is a non-parametric probability density estimate using an adaptive 
Epanechnikov kernel of the background corrected GC color distribution. The dashed red 
lines mark the bootstrapped 90% confidence limits. The solid black line represents the prob- 
ability density estimate of the uncorrected distribution. Two peaks are clearly seen at g— I 
= 0.85 and g—I — 1.35. Note that the distribution of red clusters has a broad tail. There 
is a very faint signature of the redder, metal-rich peak that is seen in normal ellipticals at 
g — I ~ 1.8. This suggests that one or both of the progenitors had only a small bulge. 
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Fig. 6. — Time evolution of the g — I color index (top panel) and absolute mag- 
nitude M v per lx 10 6 M M of sing le-burst stellar population mo dels using GA LEV 
( lAnders fc Fritze-v. Alvenslebenl 120031 1 . Model curves are plotted for a (jSalpeterlll955l ) IMF 
and metallicities as indicated in the figure. Ages of the three stellar populations are indicated 
by 0.4 Gyr, 1-1.5 Gyr and 14 Gyr for the youngest, the intermediate and the old population 
respectively. The evolution of the V magnitude, predicted from the same models is shown in 
the bottom panel and is used to compute the fading of various clusters to an age of 14 Gyr. 
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Fig. 7. — The projected spatial distribution of the clusters on the ACS image. "X" denotes 
the galaxy center. The red and blue filled triangles stand for the red and the blue clusters, 
respectively. The open red and blue squares are the red and blue objects in the background 
region where the surface number density of clusters is approximately constant. 
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Fig. 8. — J image of the inner region of AM 0139-655 showing the locations of the globular 
cluster candidates in the inner region. The dark regions denote redder colors, and the light 
regions denote bluer colors. The blue, red and magenta circles represent blue clusters, red 
clusters and clusters with g — I between 0.95 and 1.15. The bright region in the center 
is the nucleus of AM 0139-655. Note the presence of several clusters along the blue ring 
surrounding the center. 
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Fig. 9.— AM 0139-655 through the F814W filter. The frame is 84"x84". The directions on 
the sky are same as in Fig. 1. The blue and the red circles represent the blue and the red 
clusters respectively. Clusters with / brighter than 24 mag are shown by the big red circles 
and clusters with / fainter than 24 mag are shown by the small red circles. Note that the 
bright shell to the East of the galaxy is populated exclusively by blue clusters. 



-34- 




I 



Fig. 10. — The luminosity function of clusters in AM 0139-655 in the / band. The dotted 
histograms mark uncorrected (observed) LFs and the hatched histograms mark LFs corrected 
for completeness and background contamination. Completeness- weighted error bars are also 
plotted. The top panel shows the combined luminosity function of all the clusters. The 
middle panel shows the luminosity function of the blue clusters and the power law fit, giving 
an index a ~ -1.8. The lower panel shows the luminosity function of the red clusters. 
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Fig. 11. — The left panel is the luminosity function (histogram) of the red clusters. The 
curves a,b,c and d stand for estimated Gaussians representing the old metal-poor GC popu- 
lation with specific frequencies of 3.0, 1.0, 0.5 and 0.25 respectively. The histograms in the 
right panel shows the luminosity functions of the clusters obtained after subtracting the con- 
tribution from each of the gaussians in the left panel. The solid lines in the right panels are 
power-law fits to the remaining clusters. Note that these clusters can be fit by a power-law 
(solid line) with —1.6 < a < —1.85 as expected for young GCs. See Sect. 6.2 for a detailed 
discussion. 
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Fig. 12. — The CMD of the clusters associated with AM 0139-655. Superimposed are curves 
of constant mass (in units of M Q ). All curves are computed for solar metallicity with ages 
ranging from 0.3 Gyr at the leftmost point to 14 Gyr at the right extreme for each curve. 
The masses are the masses of individual clusters after fading them to 14 Gyr and which 
have current color and magnitudes as specified in the plot. Note that AM 0139-655 harbors 
several very massive red clusters. It is clear that several of these clusters will end up with 
masses equal to or greater than that of uj Cen. 



